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1 Introduction 



Although the purpose of this workshop is to discuss the structures in the universe, 
which are inhomogeneous, homogeneous models have been used in considering many 
of the cosmological issues raised in that discussion, so I have also included in this sur- 
vey the anisotropic homogeneous models and their implications. Only exact solutions 
will be covered: other speakers at Pont d'Oye (e.g. Bardeen, Brandenberger, Dunsby 
and Ellis) gave very full discussions of perturbation theory. Here I continue my pre- 
vious practice (MacCallum 1979, 1984), by using the mathematical classification of 
the solutions as an overall scheme of organization; the earlier reviews give additional 
details and references. (A survey organized by the nature of the applications is to 
appear in the proceedings of Dennis Sciama's 65th birthday meeting.) Other useful 
reviews are: Ryan and Shepley (1975) on homogeneous anisotropic models; Krasinski 
(1990) on inhomogeneous models; and Verdaguer (1985, 1992) on models of solitonic 
character. 



In section 2 I will consider the spatially-homogeneous but anisotropic models. 
These are the Bianchi models, in general, the exceptions being the Kantowski-Sachs 
models with an S 2 x R 2 topology. Such models could be significant in understanding 
the background in which structure is formed, but they do not themselves model that 
structure. However, I will include here some remarks about inhomogeneous models 
which are closely related to calculations done with Bianchi models. Then in section 
3 I will consider the inhomogeneous models, which fall into several classes. They 
can be used both as local models of structure and as possible global models of the 
background in which structure forms (and are in some cases used for both purposes 
simultaneously). A final section attempts a synthesis and makes some summarizing 
remarks. 



What is it that a cosmological model should explain? There are the following 



main features: 

[1] Lumpiness, or the clumping of matter. The evidence for this is obvious. 
[2] Expansion, shown by the Hubble law. 

[3] Evolution, shown by the radio source counts and more recently by galaxy 
counts. 

[4] A hot dense phase, to account for the cosmic microwave background radiation 
(CMWBR) and the abundances of the chemical elements. 

[5] Isotropy, shown to a high degree of approximation in various cosmological ob- 
servations, but especially in the CMWBR. 

[6] Possibly, homogeneity. (The doubt indicated here will be explained later.) 

[7] The numerical values of parameters of the universe and its laws, such as the 
baryon number density, the total density parameter Q, the entropy per baryon, 
and the coupling constants 

[8] (Perhaps) such features as the presence of life. 

Originally, the standard big-bang models were the Friedman-Lemaitre- Robertson- 
Walker (FLRW) models characterized as: 

[1] Isotropic at all points and thus necessarily. . . 

[2] Spatially-homogeneous, implying Robertson- Walker geometry. 

[3] Satisfying Einstein's field equations 

[4] At recent times (for about the last 10 10 years) pressureless and thus governed 
by the Friedman-Lemaitre dynamics. 

[5] At early times, radiation-dominated, giving the Tolman dynamics and a ther- 
mal history including the usual account of nucleogenesis and the microwave 
background. 

To this picture, which was the orthodox view from about 1965-80, the last decade 
has added the following extra orthodoxies: 

[6] 0=1. Thus there is dark matter, for which the Cold Dark Matter model was 
preferred. 

[7] Inflation - a period in the early universe where some field effectively mimics 
a large cosmological constant and so causes a period of rapid expansion long 
enough to multiply the initial length scale many times. 
[8] Non-linear clustering on galaxy cluster scales, modelled by the A^-body simula- 
tions which fit correlation functions based on observations, 
and also added, as alternatives, such concepts as cosmic strings, GUTs or TOEs[] and 
so on. 

The standard model has some clear successes: it certainly fits the Hubble law, the 
1 Why so anatomical? 



source count evolutions (in principle if not in detail), the cosmic microwave spectrum, 
the chemical abundances, the measured isotropies, and the assumption of homogene- 
ity. Perhaps its greatest success was the prediction that the number of neutrino 
species should be 3 and could not be more than 4, a prediction now fully borne out 
by the LEP data. 



However, the model still has weaknesses flMacCallum, 1987H . For example, the true 
clumping of matter on large scales, as shown by the QDOT data ||5aundcrs et ai, 1991 



and the angular correlation functions of galaxies [[Maddox et ai, 1990|| , is too strong 
for the standard cold dark matter account^. The uniformity of the Hubble flow is 
under question from the work of the "Seven Samurai" [[Lynden-Bell et ai, 1988|| and 
others. The question of the true value of Q has been re-opened, partly because theory 
has shown that inflation does not uniquely predict Q — 1 (cf. Ellis' talk at Pont 
d'Oye) and partly because observations give somewhat variant values. Some authors 
have pointed out that our knowledge of the physics valid at nucleogenesis and before 
is still somewhat uncertain, and that we should thus retain some agnosticism towards 
our account of those early times. 



Finally, we should recognize that our belief in homogeneity on a large scale has 
very poor observational support. We have data from our past light cone (and those 
of earlier human astronomers) and from geological records ||Hoyle, 1962|| . Studying 
spatial homogeneity requires us to know about conditions at great distances at the 
present time, whereas what we can observe at great distances is what happened a 
long time ago, so to test homogeneity we have to understand the evolution both of 
the universe's geometry and of its matter content^. Thus we cannot test homogeneity, 
only check that it is consistent with the data and our understanding of the theory. 
The general belief in homogeneity is indeed like the zeal of the convert, since until 
the 1950s, when Baade revised the distance scale, the accepted distances and sizes of 
galaxies were not consistent with homogeneity. 

These comments, however, are not enough to justify examination of other models. 
Why do we do that? The basic reason is to study situations where the FLRW models, 
even with linearized perturbations, may not be adequate. Three types of situation 
come to mind: the fully non-linear modelling of local processes; exploration of the 
uniqueness of features of the FLRW models; and tests of the viability of non-FLRW 
models. The uniqueness referred to here may lie in characteristics thought to be 
peculiar to the standard model; in attempted proofs that no model universe could 



2 These discoveries made it possible for disagreement with the 1980s dogmatism on such matters 
to at last be listened to. 

3 Local measures of homogeneity merely tell us that the spatial gradients of cosmic quantities are 
not too strong near us. 



be anisotropic or inhomogeneous, by proving that any strong departures from the 
standard model decay away during evolution; or in comparisons with observation, to 
show that only the standard models fit. 



Some defects of the present survey should be noted. One is that the matter content 
is generally assumed to be a perfect fluid, although this is strictly incompatible with 
the other assumed physical properties. Attempting to remedy this with some other 
mathematically convenient equation of state is not an adequate response; one must try 
to base the description of matter on a realistic model of microscopic physics or ther- 
modynamics, and few have considered such questions | Bradley and Svicstins, 1984 , 



Salvati et al., 1987, Bona and Coll, 1988, Romano and Pavon, 1992|. 



A second limitation is that we can only explore the mathematically tractable 
subsets of models^ which may be far from representative of all models. To avoid 
this restriction, we will ultimately have to turn to numerical simulations, including 
fully three-dimensional variations in the initial data. Some excellent pioneering work 
has of course been done, e.g. Anninos et al. (1991b), but capabilities are still limited 
(for example Matzner (1991) could only use a space grid of 31 3 points and 256 time 
steps). Moreover, before one can rely on numerical simulations one needs to prove 
some structural stability results to guarantee that the numerical and exact answers 
will correspond. 

As a final limitation, in giving this review I only had time to mention and discuss 
some selected papers and issues, not survey the whole vast field. For his mammoth 
survey of all inhomogeneous cosmological models which contain, as a limiting case, the 
FLRW models, Krasinski now has read about 1900 papers (as reported at the GR13 
conference in 1992)0. Thus the bibliography is at best a representative selection from 
many worthy and interesting papers, and authors whose work is unkindly omitted 
may quite reasonably feel it is unrepresentative. In particular, I have not attempted 
to cover the higher-dimensional models discussed by Demaret and others. 



2 Spatially-homogeneous anisotropic models 



4 Kramcr et al. (1980) provides a detailed survey of those classes of relativistic spacetimes where 
the Einstein field equations are sufficiently tractable to be exactly solved. 

5 The survey is not yet complete and remains to be published, but interim reports have appeared 
in some places, e.g. Krasinski (1990). 



2.1 Metrics and field equations 



As already mentioned, this class consists of the Bianchi and Kantowski-Sachs models. 
They have the advantage that the Einstein equations reduce to a system of ordinary 
differential equations, enabling the use of techniques from dynamical systems theory, 
and there is thus again a vast literature, too big to fully survey here. 



The Bianchi models can be defined as spacetimes with metrics 
ds 2 = -dt 2 + g a p(t)(e a ll dx fl )(e t3 u dx u ) 
where the corresponding basis vectors {e a } obey 

[e a , ep] = C 7 Q/3 e 7 

in which the C's are the structure constants of the relevant symmetry group. The 
different Bianchi-Behr types I-IX are then defined (see e.g. Kramer et al. (1980)) by 
algebraic classification of these sets of structure constants. 



The Kantowski-Sachs metric is 

ds 2 = -dt 2 + a 2 (t) dx 2 + b 2 (t)(d6 2 + sin 2 9 d(j) 2 ). 

(The other metric given in the original paper of Kantowski and Sachs was in fact a 
Bianchi metric, as pointed out by Ellis.) 



The adoption of methods from the theory of dynamical systems has considerably 
advanced the studies of the behaviour of Bianchi models, beginning in the early 70s 
with the discussion of phase portraits for special cases ||Collins, 1971|| . Subsequently, 



more general cases were discussed using a compactified phase space. In the last decade 
these methods have been coupled with the parametrization of the Bianchi mod- 
els using automorphism group variables | pollins and Hawking, 1973| , [Harvey, 1979| , 
Jantzen, 1979j ftiklos, 198(4 |Roque and Ellis, 1985| , |Jaklitsch, l987fl . 



The automorphism group can be briefly described as follows. Take a transforma- 



tion 



M c 



This is an automorphism of the symmetry group if the {e a } obey the same commu- 
tation relations as the {e a }. The matrices M are time-dependent and are chosen 
so that the new metric coefficients g a p take some convenient form, for example, 
become diagonal. The real dynamics is in these metric coefficients. This idea is 



present in earlier treatments which grew from Misner's methods for the Mixmaster 
Ryan and Shepley, 1975| but unfortunately the type IX case was highly mis- 



case 



leading in that for Bianchi IX (and no others except Bianchi I) the rotation group is 
an automorphism group. 

The compactification of phase space, introduced for general cases by S.P. Novikov 
and Bogoyavlenskii (see Bogoyavlenskii (1985)) entailed the normalization of config- 
uration variables to lie within some bounded region, which was then exploited by 
(a) finding Lyapunov functions, driving the system near the boundaries of the phase 
space and (b) using analyticity, together with the behaviour of critical points and 
separatrices, to derive the asymptotic behaviour. 

Three main groups have developed these treatments: Bogoyavlenskii and his col- 
leagues (op. cit); Jantzen, Rosquist and collaborators (e.g. Jantzen (1984), Rosquist 
et al. (1990)) who have coupled the automorphism variables with Hamiltonian treat- 
ments in a powerful formalism; and Wainwright and colleagues (e.g. Wainwright and 
Hsu (1989)) who have used a different, and in some respects simpler, set of automor- 
phism variables, which are well-suited for studying asymptotic behaviour because 
their limiting cases are physical evolutions of simpler models rather than singular 
behaviours. Similar ideas can be used for the Kantowski-Sachs models too. As well 
as qualitative results, some of them described below, these methods have enabled 
new exact solutions to be found, and some general statements about the occurrence 
of these solutions to be made. 



Many of the geometrical properties of Bianchi cosmologies can be carried over to 
cases where the 3-dimensional symmetry group (which is still classifiable by Bianchi 
type) acts on timelike surfaces. A number of authors have considered such metrics, 
for example Harness (1982) and myself and Siklos (1992). Although of less interest, 
since they do not evolve in time, than the usual Bianchi models, some of these models 
reappear as (spatially) inhomogeneous static or stationary models below. 

Since the present-day universe is not so anisotropic that we can readily detect its 
shear and vorticity, the Bianchi models can be relevant to cosmology only as models 
of asymptotic behaviour, in the early or late universe, or over long time scales, such 
as the time since the "last scattering" . They have also been used, in these contexts, 
as approximations in genuinely inhomogeneous universes, and one has to be careful 
to distinguish the approximate and exact uses. 



2.2 Asymptotic behaviour: the far past and future 



The earliest use of anisotropic cosmological models to study a real cosmological prob- 
lem was the investigation by Lemaitre (1933) of the occurrence of singularities in 
Bianchi type I models. The objective was to explore whether the big-bang which 
arose in FLRW models was simply a consequence of the assumed symmetry: it was 
of course found not to be. 



One can argue that classical cosmologies are irrelevant before the Planck time, 
but until a theory of quantum gravity is established and experimentally verified (if 
indeed that will ever be possible) there will be room for discussions of the behaviour 
of classical models near their singularities. 



In the late 1950s and early 60s Lifshitz and Khalatnikov and their collaborators 
showed (a) that singularities in synchronous coordinates in inhomogeneous cosmolo- 
gies were in general 'fictitious' and (b) that a special subclass gave real curvature 
singularities, with an asymptotic behaviour like that of the Kasner (vacuum Bianchi 
I) cosmology [[Lifshitz and Khalatnikov, 1963|1 - From these facts they (wrongly) in- 
ferred that general solutions did not have singularities. This contradicted the later 
singularity theorems (for which see Hawking and Ellis (1973)), a disagreement which 
led to the belief that there were errors in LK's arguments. They themselves, in col- 
laboration with Belinskii, and independently Misner, showed that Bianchi IX models 
gave a more complicated, oscillatory, behaviour than had been discussed in the ear- 
lier work, and Misner christened this the 'Mixmaster' universe after a brand of food 
mixer. The broad picture of the roles of the Kasner-like and oscillatory behaviours 
has been borne out by the more rigorous studies by the methods described in the 
previous section. There is also an interesting and as yet incompletely explored result 
that after the oscillatory phase many models approximate one of a few particular 
power-law (self-similar) solutions Pogoyavlenskii, 1985| . 



The detailed behaviour of the Mixmaster model has been the subject of still- 
continuing investigations: some authors argue that the evolution shows ergodic and 
chaotic properties, while others have pointed out that the conclusions depend cru- 
cially on the choice of time variable [ {Harrow, 1982| , [Burd et at., 1990] , |Berger, 1991 



Numerical investigations are tricky because of the required dynamic range if one is 
to study an adequately large time-interval, and the difficulties of integrating chaotic 
systems. 



The extension of these ideas to the inhomogeneous case, by Belinskii, Lifshitz 
and Khalatnikov, has been even more controversial, though prompting a smaller 



literature. It was strongly attacked by Barrow and Tipler (1979) on a number of 
technical grounds, but one can take the view that these were not as damaging to the 
case as Barrow and Tipler suggested ||Belinskii et at., 1980| , [MacCallum, 1982|| . In- 
deed the 'velocity-dominated' class whose singularities are like the Kasner cosmology 
have been more rigorously characterized and the results justified [[Eardley et a/., 1971 



Holmes et a/., 199"0|] . Sadly this does not settle the more general question, and at- 
tempts to handle the whole argument on a completely rigorous footing^ have so far 
failed. 



General results about singularity types have been proved. The 'locally extendible' 
singularities, in which the region around any geodesic encountering the singularity 
can be extended beyond the singular point, can only exist under strong restrictions 
piarke, 1976|| , while the 'whimper' singularities ||King and Ellis, 1973|| , in which cur- 
vature invariants remain bounded while curvature components in some frames blow 
up, have been shown to be non-generic and unstable [[Siklos, 1978|| . Examples of 
these special cases were found among Bianchi models, and both homogeneous and 
inhomogeneous cosmologies have been used as examples or counter-examples in the 
debate. 



A further stimulus to the study of singularities was provided by Penrose's con- 
jecture that gravitational entropy should be low at the start of the universe and this 
would correspond to a state of small or zero Weyl tensor [[Penrose, 1979|, [1'od, 199~2" 



Many authors have also considered the far future evolution (or, in closed models, 
the question of recollapse, whose necessity in Bianchi IX models lacked a rigorous 
proof until recently [[Lin and Wald, 19911 ). From various works [ MacCallum, 19 71 



|Collins and Hawking, 1973| , [Barrow and Tipler, 1978| one finds that the homogeneous 
but anisotropic models do not in general settle down to an FLRW-like behaviour 
but typically generate shears of the order of 25% of their expansion rates; see also 
Uggla et ai, 199l| . From the dynamical systems treatments, it is found that cer- 



tain exact solutions (which in general have self-similarity in time) act as attractors 
of the dynamical systems in the future ||Wainwright and Hsu, 1989| . (All such ex- 
act solutions are known: see Hsu and Wainwright (1986) and Jantzen and Rosquist 
(1986).) 



This last touches on an interesting question about our account of the evolution 
of the universe: is it structurally stable, or would small changes in the theory of 
the model parameters change the behaviour grossly? Several instances of the latter 
phenomenon, 'fragility', have recently been explored by Tavakol, in collaboration with 

6 One of them made by Smallwood and myself. 



Coley, Ellis, Farina, Van den Bergh and others 



Coley and Tavakol, 1992 



2.3 Long time effects: the cosmic microwave background 



To test the significance of the observed isotropy of the CMWBR, many people in 
the 1960s and 70s computed the angular distribution of the CMWBR temperature in 
Bianchi models (e.g. Thorne (1967), Novikov (1968), Collins and Hawking (1972), and 
Barrow et al. (1983)). These calculations allow limits to be put on small deviations 
from isotropy from observation, and also enabled, for example, the prediction of 'hot 
spots' in the CMWBR in certain Bianchi models, which could in principle be searched 
for, if there were a quadrupole component, as there is in the COBE data (though 
perhaps not for this reason), to see if the quadrupole verifies one of those models. 



Similar calculations, by fewer people, considered the polarization RRees, 1968 



Anile, 1974 , Tolman and Matzner, 1984| and spectrum [ Rees, 1968 , Rasband, 1971 



More recently still, work has been carried out on the microwave background in some 
inhomogeneous models [ 5aez and Arnau, 1990 ] . It has been shown that pure rotation 
(without shear) is not ruled out by the CMWBR fObukhov, 1992|| , but this result may 
be irrelevant to the real universe where shear is essential to non-trivial perturbations 



[Goode, 1983 


, Dunsby, 1992 


; in 


very restricted class 


Ellis, 1967 



An example of the problem with assuming a perfect fluid is that in Bianchi 
models, as soon as matter is in motion relative to the homogeneous surfaces (i.e. 
becomes 'tilted') it experiences density gradients which should lead to heat fluxes 
[Bradley and Sviestins, 1984|| : similar remarks apply to other simple models. Such 
models have recently been used to fit the observed dipole anisotropy in the CMWBR 
firmer, 1992||, though other explanations seem to me more credible. 



2.4 Early universe effects 



Galaxy formation in anisotropic models has been studied to see if by this means one 
could overcome the well-known difficulties in FLRW models (without inflation), but 
with negative results [[Perko et al., 1972]] . 



A similar investigation was to see if the helium abundance, as known in the 1960s, 
could be fitted better by anisotropic cosmologies than by FLRW models, which at the 



time appeared to give discrepancies. The reason this might happen is that anisotropy 
speeds up the evolution between the time when deuterium can first form, because 
it is no longer dissociated by the photons, and the time when neutrons and protons 
are sufficiently sparse that they no longer find each other to combine. Hawking and 
Tayler (1966) were pioneers in this effort, which continued into the 1980s but suffered 
some mutations in its intention. 



First the argument was reversed, and the good agreement of FLRW predictions 
with data was used to limit the anisotropy during the nucleogenesis period (see e.g. 
Barrow (1976), Olson (1977)). Later still these limits were relaxed as a result of con- 
sidering the effects of anisotropic neutrino distribution functions [[Rothman and Matzner, 1982|| 
and other effects on reaction rates llJuszkiewicz et ai, 198311 . It has even been shown 



[[Matravers et ai, 1984[ Parrow, 1984|| that strongly anisotropic Bianchi models, not 
obeying the limits deduced from perturbed FLRW models, can produce correct ele- 
ment abundances, though they may violate other constraints [[Matravers and Madsen, 1985| , 
Matravers et a/., 1985 



3 Inhomogeneous cosmologies 



3.1 Self-similar models 



Some of the self-similar models, especially those relevant to modelling structure for- 
mation, are reviewed in much greater detail in a complementary talk by Carr, so I 
will give here only a few details of other cases. 

The geometry of the self-similar models first considered in cosmology is somewhat 
like that of the Bianchi models, except that one of the isometries is replaced by a 
homothety, that is to say by a vector field satisfying 



f ( ;b) = 2kg. 



a b 



where A; is a constant. This class, where the homothety and two independent sym- 
metries act, was considered by a number of authors fEardley, 1974J , |Luminet, 1978 



Wu, 198IT , [Hanquin and Demaret, 1984I H, and many details, parallel in nature to 



those covered by the detailed studies of Bianchi models, can be found in those works. 



7 Due to Western confusion over Chinese name order, Wu Zhong-Chao is sometimes incorrectly 
referred to as W.Z. Chao rather than Wu, Z-C. 



More recently Wainwright, Hewitt and colleagues [ [Hewitt et al., 1955 , [Hewitt and Wainwright, 
Hewitt et at., 199~T| have considered cases where the homothety has a timelike rather 



than spacelike generator. Like the former class, these solutions are in fact special cases 
of "G2 solutions" (discussed below) with perfect fluid matter content. It is found that 
the spatial variations can be periodic or monotone; the asymptotic behaviour may 
be a vacuum or spatially homogeneous model; the periodic cases are unstable to 
increases in the anisotropy; and the singularities can be acceleration-dominated. 



3.2 Spherically symmetric models 



These have a metric 

d s 2 = -e 2v dt 2 + e 2X dr 2 + R 2 (d6 2 + sin 2 9 d<j> 2 ) 

where v, A and R are functions of r and t. The precise functional forms in the metric 
depend on the choice of coordinates and the additional restrictions assumed. It should 
be noted that there are so few undetermined functions that a sufficiently-complicated 
energy-momentum will fit a totally arbitrary choice of the remaining functions: in my 
view this should not be regarded as a solution, since no equation is actually solved! 



Some important subcases have been studied, notably: 
[1] The dust (pressureless perfect fluid) cases, originally studied by Lemaitre, but 

usually named after Tolman and Bondi; 
[2] McVittie's 1933 solution representing a black hole in an FLRW universe; 
[3] The "Swiss cheese" model constructed by matching a Schwarzschild vacuum 

solution inside some sphere to an exterior FLRW universe; 

Shearfree fluid solutions ||Wyman, 1946| , |Kustaanheimo and Qvist, 1948| , [Stephani, 1983| . 
McVittieT 



[4] 

[5] Self-similar solutions, discussed in Carr's contribution at Pont d'Oye. 



T981 



Spherically symmetric models, especially Tolman-Bondi, have often been used to 
model galactic scale inhomogeneities, in various contexts. Galaxy formation has been 
studied (e.g. Tolman (1934), Carr and Yahil (1990)): Meszaros (1991) developed a 
variation on the usual approach by considering the shell-crossings, with the aim of 
producing "Great Wall" like structures, rather than the collapse to the centre pro- 
ducing a spherical cluster or galaxy. Some authors have used spherically symmetric 
lumps to estimate departures from the simple theory of the magnitude-redshift re- 
lations based on a smoothed out modelQ (e.g. Dyer (1976), Kantowski (1969a) and 

8 The point is that the beams of light we observe are focussed only by the matter actually inside 
the beam, not the matter that would be there in a completely uniform model. 



Newman (1979)): note that these works show that the corrections depend on the 
choice of modelling, since Newman's results from a McVittie model differ from the 
ones based on Swiss cheese models. The metrics also give the simplest models of 
gravitational lenses^ and have also been used to model the formation of primordial 
black holes |Carr and Hawking, 1974 |. 



On a larger scale, inhomogeneous spherical spacetimes have been used to model 
clusters of galaxies [ Kantowski, 1969b | , variations in the Hubble flow due to the super- 



cluster Mavrides, 1977 



the evolution of cosmic voids fSato, 1984] , |Hausman et al, 



1983, 



[Bomior and Chamorro, 1990 and 1991 1 , the observed distribution of galaxies and sim- 



ple hierarchical models of the universe fBonnor, 1972| , [Wesson, 1978| , |Wesson, 1979 



[Ribciro, 1992a . Most of this work used Tolman-Bondi models, sometimes with dis- 
continuous density distributions. 



Recent work by Ribeiro (1992b), in the course of an attempt to make simple 
models of fractal cosmologies using Tolman-Bondi metrics, has reminded us of the 
need to compare data with relativistic models not Newtonian approximations. Taking 
the Einstein-de Sitter model, and integrating down the geodesies, he plotted the 
number counts against luminosity distances. At small distances, where a simple 
interpretation would say the result looks like a uniform density, the graph is irrelevant 
because the distances are inside the region where the QDOT survey shows things are 
lumpy [[Saunders et at., 1991 1, while at greater redshifts the universe ceases to have a 
simple power-law relation of density and distance. Thus even Einstein-de Sitter does 
not look homogeneous! 



One must therefore ask in general "do homogeneous models look homogeneous?" . 
Of course, they will if the data is handled with appropriate relativistic corrections, but 
to achieve such comparisons in general requires the integration of the null geodesic 
equations in each cosmological model considered, and, as those who have tried it 
know, even when solving the field equations is simple, solving the geodesic equations 
may not be. 



Many other papers have considered spherically symmetric models, but there is 
not enough space here to review them all, so I will end by mentioning a jeu d 'esprit 
in which it was shown that in a "Swiss cheese" model, made by joining two FLRW 
exteriors at the two sides of a Kruskal diagram for the Schwarzschild solution, one 
can have two universes each of which can receive (but not answer) a signal from the 
other ||Sussman, 1985|| . 

9 The very detailed modern work interpreting real lenses to study various properties of individual 
sources and the cosmos mostly uses linearized approximations. 



3.3 Cylindrically symmetric and plane symmetric (static) 
models 



These have been used to model cosmic strings and domain walls. One should note 
that locally the metrics may be the same for these two cases, the difference lying in 
whether there is or is not a Killing vector whose integral curves have spatial topology 
S 1 . Plane symmetric metrics should have a rotational symmetry in the plane, but 
to add to the possible confusions some authors use the term "plane" for solutions 
without such a rotation: the term "planar" would be a useful alternative. 

The usual (though not the only) form for the cylindrically symmetric metrics is 

ds 2 = -T 2 dt 2 + R 2 dr 2 + Z 2 dz 2 + 2Wdzd(j) + <5> 2 d<p 2 

where T, R, Z, W and $ depend on r (and, in the non-stationary case, t) and (ft is 
periodic, and, for the plane symmetric case, 

ds 2 = -T 2 dt 2 + R 2 dr 2 + X 2 (dx 2 + x 2 d(f) 2 ) 

where T, R and X are functions of r (and perhaps t). The static cases all belong in 
Harness's (1982) general class. 

Plane symmetric models, usually static, solutions have been used to model domain 
walls ||Vilenkin, 19831 , [Ipser and Sikivie, 1984| , |Goetz, 19901 |Wang, 1991b| H- The cylin- 



drically symmetric models have been used for cosmic strings, starting with the work 
of Gott, Hiscock and Linet in 1985. These studies have usually been done with static 
strings^, and have considered such questions as the effects on classical and quantum 
fields in the neighbourhood of the string. 



3.4 G2 cosmologies 



I use the above title as a general name for all cosmological metrics with two spacelike 
Killing vectors (and hence two essential variables). The cylindrical and plane metrics, 
and many of the Bianchi metrics, are special cases of G2 cosmologies. 

G2 cosmologies admit a number of specializations, such as: 

10 Note that since the sources usually have a boost symmetry in the timelike surface giving the wall, 
corresponding solutions have timelike surfaces admitting the (2+l)-dimcnsional de Sitter group. 

11 There is some controversy about whether these can correctly represent strings embedded in an 
expanding universe flClarke et al, 1990| . 



[1] the Killing vectors commute; 

[2] the orbits of the G<i are orthogonal to another set of 2- dimensional surfaces V2] 
[3] the Killing vectors individually are hypersurface-orthogonal; 
[4] the matter content satisfies conditions allowing generating techniques. 
Among the classes of metrics covered here are colliding wave models, cosmologies with 
superposed solitonic waves, and what I call "corrugated" cosmologies with spatial 
irregularities dependent on only one variable. 

The metrics where the Killing vectors do not commute have been very little stud- 
ied: it is known they cannot admit orthogonal V2 if the fluid flows orthogonal to the 
group surfaces (unless they have an extra symmetry) and that if the fluid is thus 
orthogonal it is non-rotating ||Bugallio, 1987| . |van den Bcrgh, 1988 ] . So we now take 



only cases where the Killing vectors commute. 

The case without orthogonal V2 has also been comparatively little studied, but 
recently some exact solutions which have one hypersurface-orthogonal Killing vector 
and in which the metric coefficients are separable, have been derived and studied 
van den Bcrgh et at., 1991 , van den Bcrgh, 1991|| . One class consists of metrics of 



the form 

ds 2 = e 2{ ^ K+k \-dt 2 + dx 2 ) + e 2 ^ s+s \e F+i dy) 2 + {e^ F+ ^6) 2 } 

where: K, S and F depend on t; k, s and / depend on x; 8 = dz + 2uj dx; and uj 
depends on t and x. Some perfect fluid solutions are known explicitly but usually 
turned out to be self-similar, with big-bang singularities of the usual types. The "stiff 
fluid" (7 = 2) is a special case, discussed in detail by van den Bergh (1991). Most 
of the solutions have singularities at finite spatial distances or can be regarded as 
inhomogeneous perturbations of the Bianchi VI-\ models. 

The cases with orthogonal V2 were classified by Wainwright (1979;1981), and a 
number of specific examples are known (e.g. Wainwright and Goode (1980); Kramer 
(1984)). A recent solution found by Senovilla (1990) attracted much attention, be- 
cause it is non-singular [phinea et ai, 1991|| , evading the focussing conditions in the 



singularity theorems by containing matter that is too diffuse: it is closely related 
to an earlier solution of Feinstein and Senovilla (1989)Q. The metrics investigated 
in this class generally have Kasner-like behaviour near the singularity (though some 
have a plane- wave asymptotic behaviour ||Wainwright, 1983| ) and become self-similar 



or spatially homogeneous in the far future. 

Finally we come to the most-studied class, those where the generating techniques 



12 Some recent work has g iven generalizations of these solutions | Ruiz and Senovilla, 1992[ 
van den Bcrgh and Skca, 1992f] ; also S.W. Goode at GR13 (unpublished). 



are applicable. The matter content must have characteristic propagation speeds equal 
to the speed of light, so attention is restricted to vacuum, electromagnetic, neu- 
trino and "stiff fluid" (or equivalently, massless scalar field with a timelike gradient) 
cases. However, FLRW fluid solutions can be obtained by using the same methods in 
higher- dimensions and using dimensional reduction. There are useful reviews cover- 
ing the cosmological, cylindrical, and colliding wave sub-classes [ Carmeli et al, 1981 



Verdaguer, 1985| , |Verdaguer, 1992| , |Ferrari, 1990| , [Griffiths, 199 1|] . The metrics can be 



written in a form covering also the related stationary axisymmetric metrics as 

ds 2 = ef AB dx A dx B + 5e 27 ((dx 4 ) 2 - e(dx 3 ) 2 )/f 
where A, B take values 1, 2 and the values of Jab can be written as a matrix 

{ -fu fu 2 + e(x 3 ) 2 // ) 

The case 5 = — e = 1 gives the stationary axisymmetric metrics, the case 5 = e = 1 
the cylindrical cases and e = —5 = 1 the cosmological cases. Physically these classes 
differ in the timelike or spacelike nature of the surfaces of symmetry and the nature 
of the gradient of the determinant of the metric in those surfaces. 

Some studies have focussed on the mathematics, showing how known vacuum so- 
lutions can be related by solution-generating techniques [[Kitchingham, 19"84| , while 



others have concentrated on the physics of the evolution and interpretative issues. 
The generating techniques use one or more of a battery of related methods: Backlund 
transformation, inverse scattering, soliton solutions and so on. One interesting ques- 
tion that has arisen from recent work is whether solitons in relativity do or do not 
exhibit non-linear interactions: Boyd et al. (1991), in investigations of solitons in 
a Bianchi I background, found no non-linearity, while Belinskii (1991) has claimed 
there is a non-linear effect (see also Verdaguer (1992)). 

The applications in cosmology, which have generated far too many papers to list 
them all here, have been pursued by a number of groups, notably by Carmeli, Charach 
and Feinstein, by Verdaguer and colleagues, by Gleiser, Pullin and colleagues, and 
Belinski, Curir and Francaviglia, with important contributions by Ibanez, Kitching- 
ham, Yurtsever, Ferrari, Chandrasekhar and Xanthopoulos, Letelier, Tsoubelis and 
Wang and many others. 

One use of these metrics is to provide models for universes with gravitational 
waves. It emerges that the models studied are typically Kasner-like near the sin- 
gularity (agreeing with the LK arguments), and settle down to self-similar or spa- 
tially homogeneous models with superposed high-frequency gravitational waves at 



late times IIAdams et a/., 19821, 



Carmeli and Feinstein, 1984|, |Feinstein, 1988 



An- 



other use is to model straight cosmic strings in interaction with gravitational or 
other waves (e.g. Economou and Tsoubelis (1988), Verdaguer (1992)). One can 
also examine the gravitational analogue of Faraday rotation [[Piran and Safler, 1985|, 



|Tomimatsu, 19891 , |Wang, 1991a|1 and there are even solutions whose exact behaviour 
agrees precisely with the linearized perturbation calculations for FLRW universes 



Carmeli et al, 1983 



3.5 Other models 



Solutions with less symmetry than those above have been little explored. Following 
Krasinski one can divide the cases considered into a number of classes (in which I 
only mention a few important special subcases). 



The Szekeres-Szafron family (also independently found by Tomimura). These 
have in general no symmetries, and contain an irrotational non-accelerating 
fluid. Tolman-Bondi universes are included in this class, as are the Kantowski- 
Sachs metrics; some generalizations are known, such as the rotating inhomoge- 
neous model due to Stephani (1987). Like the G2 solutions mentioned earlier, 
some Szekeres models obey exactly the linearized perturbation equations for 



the FLRW models [|Goode and Wainwright, 1982 



Shearfree irrotational metrics [Parnes, 1973] which include the conformally flat 



fluids (Stephani 1967a, 1967b) and McVittie's spherically symmetric metric. 
Bona and Coll (1988) have recently argued that the Stephani cases can only 
have acceptable thermodynamics if the metrics admit three Killing vectors. 

3. The Vaidya-Patel-Koppar family, which represent an FLRW model contaning a 
"Kerr" solution using null radiation and an electromagnetic field. The physical 
significance of these metrics is dubious. 

4. Some other special cases such as Oleson's Petrov type N fluid solutions. 



4 Syntheses and conclusions: what have we learnt? 



Here I collect up the outcome of the work surveyed above, without repeating all 
details, and review some relevant extra references, but many interesting aspects are 



still omitted. For example, the literature covers such issues as models for interactions 
between different forms of matter, and generation of gravitational radiation. 



4.1 The classical singularity 



The occurrence of a "big-bang" in FLRW models is not just a consequence of the 
high symmetry. Its nature in general models is probably a curvature singularity, and 
the best guess so far is that the asymptotic behaviour would be oscillatory but other 
possibilities exist. The Penrose conjecture, which would be a selection principle on 
models, has been particularly developed, using exact solutions as examples, by Wain- 
wright and Goode, who have given a precise definition of the notion of an 'isotropic 
singularity' [Coode et ag 1992| , [Tod, 1992|| . 



4.2 Occurrence of inflation 



In "old" inflation in Bianchi I models, inflation need not occur ||Barrow and Turner, 1981 



but in "new" inflation it was predicted ||Stcigman and Turner, 1983|| . In a large class 
of chaotic inflation models it is also expected [ |Moss and Sahni, 1986[ |. Further pa- 
pers by a number of authors have suggested that inflation need not always occur (see 
Rothman and Ellis (1986) for some criticisms of earlier papers). 



4.3 Removal of anisotropy and inhomogeneity 

Three means of smoothing the universe have been explored over the years: the use of 
viscosity in the early universe; the removal of horizons in the Mixmaster universes; 
and removal during inflation. The first two of these ingenious suggestions are due to 
Misner. 

Attempts to smooth out anisotropics or inhomogeneities by any process obeying 
deterministic sets of differential equations satisfying Lipschitz-type conditions are 
doomed to fail, as was first pointed out by Collins and Stewart (1971) in the context 
of viscous mechanisms. The argument is simply that one can impose any desired 
amount of anisotropy or inhomogeneity now and evolve the system backwards in 
time to reach initial conditions at some earlier time whose evolution produces the 
chosen present-day values. 



The same argument also holds for inflationary models. Inflation in itself, without 
the use of singular equations or otherwise indeterminate evolutions, cannot wholly ex- 
plain present isotropy or homogeneity, although it may reduce deviations by large fac- 
tors | [Sirousse-Zia, 198"!^ , |Wald, 1983| , [Moss and Salmi, 1986| , [Kutamase et at., 1989 



Although one can argue that anisotropy tends to prolong inflation, this does not 
remove the difficulty 

Since 1981 I have been arguing a heretical view about one of the grounds for 
inflation, namely the 'flatness problem', on the grounds that the formulation of this 
problem makes an implicit and unjustified assumption that the a priori probabilities 
of values of Q is spread over some range sufficient to make the observed closeness to 
1 implausible. Unless one can justify the a priori distribution, there is no implausi- 
bilityQ |[Ellis, 1991|| . 



However, if one accepts there is a flatness problem, then there is also an isotropy 
problem, since at least for some probability distributions on the inhomogeneity and 
anisotropy the models would not match observation. Protagonists of inflation cannot 
have it both ways. Perhaps, if one does not want to just say "well, that's how the 
universe was born", one has to explain the observed smoothness by appeal to the 
'speculative era', as Salam (1990) called it, i.e. by appeal to one's favourite theory of 
quantum gravity. 

If inflation works well at early times, then inflation actually enhances the chance 
of an anisotropic model fitting the data, and since the property of anisotropy cannot 
be totally destroyed in general (because it can be coded into geometric invariants 
which cannot become zero by any classical evolution) the anisotropy could reassert 
itself in the future! (This of course will not happen if a non-zero A term persists, as 
the "cosmic no-hair" theorems show ||Wald, 1983| , [Morrow- Jones and Witt, 1988|| .) 



The Mixmaster horizon removal suggestion was shown to fail when more detailed 
computations than Misner's were made [[MacCallum, 1971], |Doroshkevich et ai, 1971 



|Chitre, 1972|| . Incidentally, one may note that inflation does not solve the original 
form of the 'horizon problem', which was to account completely for the similarity of 
points on the last scattering surface governed by different subsets of the inital data 
surface. Inflation leads to a large overlap between these initial data subsets, but not 
to their exact coincidence. Thus one still has to assume that the non-overlap regions 
are not too different. While this may give a more plausible model, it does not remove 



13 One can however argue that only fi = 1 is plausible, on the grounds that otherwise the 
quantum theory before the Planck time would have to fix a length-scale parameter much larger 
than any quantum scale, only the 0=1 case being scale-free. I am indebted to Gary Gibbons for 
this remark. 



the need for assumptions on the initial data. 



4.4 The exit from inflation 

A further interesting application of non-standard models has come in a recent attempt 
to answer the question posed by Ellis and Rothman (unpublished) of how the universe 
can choose a uniform reference frame at the exit from inflation when a truly de Sitter 
model has no preferred time axis. Anninos et al. (1991a) have shown by taking an 
inflating Bianchi V model that the answer is that the memory is retained and the 
universe is never really de Sitter. 



4.5 The helium abundance 

This is still used to set limits on anisotropy during the nucleosynthesis phase. 



4.6 The cosmic microwave background 

Observations limit the integrated effect since "last scattering" : note this can in princi- 
ple permit large but compensating excursions from FLRW. One intriguing possibility 
raised by Ellis et al. (1978) is that the observed sphere on the last scattering surface 
could lie on a timelike (hyper) cylinder of homogeneity in a static spherically sym- 
metric model. This makes the CMWBR isotropic at all points not only at the centre, 
and although it cannot fit all the other data, the model shows how careful one must 
be, in drawing conclusions about the geometry of the universe from observations, not 
to assume the result one wishes to prove. 

There is a theorem by Ehlers, Geren and Sachs (1968) showing that if a congruence 
of geodesically-moving observers all observe an isotropic distribution of collisionless 
gas the metric must be Robert son- Walker. Treciokas and Ellis (1971) have inves- 
tigated the related problem with collisions. Recently Ferrando et al. (1992) have 
investigated inhomogeneous models where an isotropic gas distribution is possible. 
These studies throw into focus a conjecture which is usually assumed, namely that an 
approximately isotropic gas distribution, at all points, would imply an approximately 
Robert son- Walker metric. (It is this assumption which underlies the arguments nor- 
mally used in analysing data like that from COBE to get detailed information on 



allowed FLRW perturbations.) 



4.7 The far future 

Anisotropy will in general become apparent, if it is present and if the cosmological 
constant A is zero: isotropy is not stable. Inhomogeneities may become significant 
even faster. 



4.8 The origin of structure 

None of the work discussed above accounts for the origin of structure, although it 
offers suitable descriptions for the evolution, or the background spacetime in which 
the evolution takes place. I feel it does, however, indicate strongly that the true origin 
lies in the perhaps unknowable situation in the Speculative Era, and the resulting 
initial conditions for the later evolution. 



4.9 A genuinely anisotropic and inhomogeneous universe? 

While I do not think one can give a definitive answer to this question, I would 
personally be very surprised if anisotropic but homogeneous models turned out to 
be anything more than useful examples. However, the status of fully inhomogeneous 
models is less clear. 

One argument is that while the standard models may be good approximations 
at present, they are unstable to perturbations both in the past and the future. The 
possible alternative pasts are quite varied, as shown above, even without considering 
quantum gravity. Similarly, the universe may not be isotropic in the far future. 
Moreover, we have no knowledge of conditions outside our past null cone, where 
some inflationary scenarios would predict bubbles of differing FLRW universes, and 
perhaps domain walls and so on. 

If the universe were FLRW, or very close to that, this means it is in a region, in the 
space of all possible models, which almost any reasonable measure is likely to say has 
very low probability (though note the earlier remarks on assignments of probabilities). 
One can only evaluate, and perhaps explain, this feature by considering non-FLRW 



models. It is noteworthy that many of the "problems" inflation claims to tackle are 
not problems if the universe simply is always FLRW. Hence, as already argued above, 
one has a deep problem in explaining why the universe is in the unlikely FLRW state 
if one accepts the arguments about probabilities current in work on inflation. 

Suppose we speculated that the real universe is significantly inhomogeneous at 
the present epoch (at a level beyond that arising from perturbations in FLRW). 
What would the objections be? There are only two relevant pieces of data, as far 
as I can see. One is the deep galaxy counts made by the automatic plate measuring 
machines, which are claimed to restrict variations to a few percent, and the other is 
the isotropy of the CMWBR. Although the latter is a good test for large lumps in 
a basically FLRW universe, one has to question (recalling the results of Ellis et al. 
(1978)) whether it really implies homogeneity. 
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